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1 Introduction 



There has been much recent work which suggests that the baryon asymmetry of 
the universe might be produced at the electroweak phase transition Q. If the 
equations of motion of electroweak theory are quantized, then anomalies arise which 
are responsible for the Chern-Simons -or baryon number changing- currents 0. 
The rate of baryon number violation is rather small for low temperatures but it 
becomes large at a temperature scale on the order of the electroweak scale of about 
100 GeV H - pi. Since baryon number is violated, a net baryon number will be 
generated if there is a mechanism that biases the rate at the electroweak phase 
transition in such a way, that as the baryon number violation shuts off in the low 
temperature phase, an asymmetry remains frozen into the system. This can happen 
in practice if the electroweak phase transition is of first order, and if there is sufficient 
CP violation at the temperature of the phase transition. 

It is expected that the electroweak phase transition is of first order § -||. In a 
first order phase transition, the conversion from one phase to another occurs through 
nucleation of the true phase in the false phase. This happens when the system is 
either supercooled or superheated. The bubbles of the true phase expand rapidly 
eating up the region of the false phase. For the electroweak phase transition, this 



true phase eventually fills the entire volume with no intermediate mixed phase [1C]. 
At the bubble surface, there is a thin wall of microscopic dimensions which separates 
the phases. It is at this bubble wall that matter is strongly out of equilibrium, and 



here the baryon asymmetry is generated [11| - [jig ]. 

To quantitatively understand the generation of the baryon asymmetry in the 
bubble wall, the effect of the bubble wall on the propagation of fermions should 
be understood. Fermions passing through the bubble or domain wall acquire mass, 
generated by the Yukawa coupling, which is proportional to the finite temperature 
vaccuum expectation value (vev) of the Higgs field. This vev is determined from 
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the equations of motion of the finite temperature effective action of the bubble. 

To simplify the problem, we will work in the approximation where the energy 
density of the two phases are degenerate. In practice this is a good approxima- 
tion, since in cosmology, the universe is expanding so slowly that the nucleation 
always begins at a temperature where the amount of supercooling is very small, 
and the energies of the two phases are degenerate. In this approximation, there is 
a one dimensional kink solution which separates the phases, and the bubble wall 
is a domain wall which separates the domains of different energy. This kink wall 
can propagate at any velocity. In practice, the wall velocity is determined by a 
complicated analysis which involves computing the effects of dissipative processes, 
and is in the range 0.1 - 0.9 of the speed of light ||, [JOJ - [21|. In our analysis, we 
will consider fermion propagation in the presence of the wall at rest. The case for 
a moving wall can be determined by Lorentz boosting to the moving wall frame. 

A plot of the domain wall is shown in Fig. 1. At x = — oo, the system is in 
the symmetric phase, that is outside the bubble. At x = +oo, the system is in 
the symmetry broken phase, that is, inside the bubble. The approximation of the 
bubble as a planar interface should be valid for bubbles which are large compared to 
a microscopic size scale. This is true for most of the evolution of bubbles produced 



in the electroweak phase transition [10]. 

The outline of this paper is as follows: in Section 2 we describe how to obtain 
the equation of motion for fermions in the presence of a domain wall in the minimal 
standard model. We obtain the Dirac equation with an effective mass proportional 
to the vev of the Higgs field. We solve this equation in Section 3, and compute 
the fermion wave function we obtain transmission and reflection coefficients. In 
Section 4, we present the normalized solutions of the Dirac equation. We finish by 
summarizing our results in section 5. 
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2 The Motion of Fermions in the Presence of Domain 
Walls 



In this section we show how to obtain simplified classical equations of motion for 
fermions in the standard model scattering from and interacting with the electroweak 
domain walls. We use the classical mean-field approximation in which the bosonic 
field operators are replaced by their classical expectation values. 

We will briefly review the notation which we will use throughout the paper. 
The Lagrangian in the standard model is invariant under SU(2) x U(l) transfor- 
mations and the fields are therefore eigenstates of weak isospin and hypercharge Y. 
The SU(2) x U(l) invariant vacuum state of the Lagrangian exists only in a high 
temperature phase above the electroweak phase transition. The symmetry is broken 
spontaneously once we cool the system below the transition temperature T c . 

The Lagrangian is 

^gaugefield ^Higgs ~i~ ^-fermion ^--Yukawa ^gf- (-Q 

The gauge field, fermion, scalar kinetic energy and gauge fixing terms are all 



written in the standard way [24]. For the potential energy of the Higgs field, we 
take 

V($) = \{&<S> -v 2 /2) 2 (2) 
so that the vacuum expectation value of the Higgs field is 

To describe kink solutions and a first order phase transition, one has to include 
the modifications of the effective potential for the Higgs field due to interactions 
with the heat bath. Including effects of both one-loop diagrams and the sum of ring 
diagrams, the effective potential can be written to a good approximation as 

V eff (<f>,T) = |(T 2 - T C V - 5T<P 3 + ^ 4 (4) 



where (j> = \/2($^ < I>) 1 / 2 . The dimensionless parameters 7 and 5 are approximately 
given by 

5 2 
1 



7 ~ — J 
' 16 s 



when Mw ~ -^z ~ Mt, where these are the masses of the W-bosons, the Z-boson 
and the top respectively. 5 is small and is responsible for generating the first order 
phase transition. At some value of T, there are two degenerate minima of this 
potential and a kink solution exists. 

Of course the above evaluation of the effective potential may not be such a 
good approximation for realistic values of the Higgs mass. If the effective potential 
is greatly modified then the kink solution, estimates for bubble wall velocity and the 
considerations we present will all have to be modified. Nevertheless, in the words 
of an old U. S. television show, it is better to light one little candle than to curse 
the darkness. 

As the scalar field evolves through the kink solution, no current is generated 
for the vector field, so that no expectation value of the vector fields is generated. 
We are therefore justified in truncating the system to only the Higgs field degrees of 
freedom. The equation of motion is just the classical Higgs field equation of motion 
with zero background vector field. 



To get the effective potential for the Higgs field into more transparent form [IS], 
we introduce the dimensionless temperature £ = (1 — (t') ) anc ^ the dimensionless 
field strength g = -^<f) and obtain 

V ef f(9) = 6T( 5 ^ ffy-g^-/) (6) 

The potential V e ff{g) is plotted in Figure 2 for different values of £. For large 
positive values of £ we are in the high temperature phase with the vev of the Higgs 
field being zero. Decreasing £ we develop at ( = 2.25 a second relative minimum 



which for Q = 2 becomes degenerate with the first one. For smaller values of (, 
the high temperature phase at g = becomes unstable with respect to the new 
absolute minimum at g ^ , the first order phase transition sets in and bubbles of 
the broken phase start to nucleate. If we supercool to £ = the system spinodally 
decomposes. 

Since the expansion rate of the universe (« 10~ n s) is rather slow compared to 
the electroweak timescale (« 10~ 26 s) we cannot strongly supercool and the phase 
transition will roughly proceed at ( ~ 2 where the minima of the effective potential 
are degenerate. The steady state solution describing the domain wall is thus given 
approximately by the solution of the Higgs field g for a transition between the 
degenerate minima at temperature £ = 2. 

The equation of motion for the kink is mathematically identical to treating the 
amplitude of the Higgs field g as a spatial coordinate for the inverted potential of 
the previous equation and taking the spatial coordinate of the kink as the time 
variable. The mechanical analog to this system is a frictionless ball rolling from the 
top of one hill through a valley to the top of another hill. For degenerate minima 
the tops of the two hills have the same height and there exists a solution where the 
ball starts at the top of one hill from rest and ends on the top of the other hill again 
at rest. Such a solution is the kink. Energy conservation for the kink reads 

E = \(f r ) 2 -V eff (g(r)) (7) 

which for a transition between degenerate minima (E = 0) and a kink bubble with 
dimensionless position x = ^j^Tr can be rewritten as 

§--,(, -2). (B) 

We can integrate the above equation implementing the boundary condition that 
at the center of the nucleated bubble (x — > +oo) the system is in the broken phase 
with a finite vev while outside the domain wall (x — > — oo) it approaches a zero 



vev in the high temperature phase, treating the domain wall as effectively flat and 
separating two semi-infinite domains. We obtain 



g = 1 + tanh(x). 



0) 



This solution, commonly called the 'kink', is depicted in figure 1. It represents the 
classical, steady state, finite temperature solution to the equations of motion. 

Given a background scalar field, we can solve the Dirac equation in the presence 
of this field. Rewriting the coordinates in terms of ordinary dimensionful length 
scales, the Dirac equation is 



where ^ is any fermion field, and £ is twice the ratio of quark to Higgs mass at 
zero temperature. The solution to this equation describes the motion of a fermion 
in the presence of a kink. To this order of approximation for the standard model, 
CP violation plays no essential role in this equation. 



citly showing how to solve this equation analytically for the fermion wave function. 
The second part is devoted to deriving the transmission and reflection coefficients 
from the wave function. 

3.1 The wave function 

We will solve the Dirac equation by looking for eigenstates of momentum in the 
plane of the wall, pt, that is the momentum perpendicular to a normal vector on 
the surface of the wall. The Dirac equation (|i~0|) is then a differential equation in one 
dimension and the only relevant variable (r) is the one along the direction normal to 




(10) 



3 Solving the Dirac equation 



In this section we study the solutions to the Dirac equation ([n]) . We start by expli- 
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the domain wall. The degrees of freedom in the plane of the domain wall effectively 
decouple. This can be used by substituting the ansatz 

5T 



¥(r, x t ,t) 

into equation (|1C|) , finding 

5T 



■ + 



^g(r))e ±i{ptXt - Et) ^(r) 



5T 



(11) 



(12) 



The + and — signs in equation ( |1 1[ ) correspond to the positive and negative energy 
solutions respectively. 

The commutator in equation (|12|) picks only a contribution from the direction 
normal to the wall and therefore it is proportional to just one of the 7 -matrices 
and we choose this to be 73. We will use the ordinary representation for gamma 
matrices where 

/ 
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V 



a 1 
-a 1 



(13) 



The eignestates of 73 are 



73 u s ± = ±iu s > 



(14) 



where s = 1, 2 and are 



/ 1 \ 



±i 
V J 



and u 4 



/ \ 
1 


V Ti ) 



(15) 



Recalling that the dimensionless radius x is x = -J=Tr, we obtain from (|l~2|), 
in dimensionless variables 



d 2 d 
{^2 + ^^(x) ~ eV(^) + e 2 }<S>(x) = 



(16) 



The dimensionless energy parameter e is given in terms of the energy E and the 
momentum parallel to the bubble wall pt by 



<2\/8T)^E*-pl (17) 
We write the spinor <£(x) in the basis ( |l5[ ) 

$(*) = £</4(*)u± (is) 

The spinors ( |l~5| ) are linearly independent and thus, writing (18) into equation ( |i~6| ) 
we get two differential equations for the functions 4>±, 

0. (19) 

To further simplify equations (|l9[), we factorize the singularities. With the 
change of variable z = i(l — tanh(x)) we substitute 

Mx)=z a (l-zf X± (z) (20) 

and examine the behavior of the resulting differential equations near the singular 
points z = and z = 1. Assuming that the functions x± are slowly varying near 
the singularities, this yields two algebraic conditions from which we get 



ax z ax 



a 



i 

+ 2 £ 



4£ 2 
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(the choice of signs for a and (3 is discussed later). With the values of a and (3 at 
hand, the differential equation for the functions x± is the hypergeometric differential 
equation 

d 2 d 
{z{\ - z )— + { c -{l + a T + b T )z)— - a T b T } X ± = (22) 

with the parameters a=p , 6=p and c given by 



C 



a + f3 + -- 

a + (3+ - + 
2a + 1. 



1 



(23) 



Each of the equations (22) has two independent solutions |2^] 

X ± a \z) = 2 F 1 (a T ,b T ,c;z) 

X± a \z) = z 1 - 2 F 1 (a T + l-c,b T + l-c,2-c;z) (24) 

The hypergeometric equations are expressed here as expansions around z = 0. The 
superscripts (±a) are explained below. 



The general solutions to equations (19) are 



(A s )^ a U { ± a) + (A s ) { + a U ( ± a) (25) 



with 



= z a (l-zf 2 F 1 (a T ,b T ,c;z) 
4 +Q) = z~ a {l-zf 2 F 1 (a T + l-c,b T + l-c,2-c;z) (26) 

where we have used (p0[). The superscripts (±a) indicate the behavior of these 
solutions at x — * +oo or correspondingly at z = since 

z ±a {l-zf e^ 2ax (27) 

z a {l-z) ±f3 e ±2f3x . (28) 

and 2 Fi(a, b, c; 0) = 1. The behavior of the solution in these limits will be needed 
in the next section to implement boundary conditions. 

We now proceed to find the fermion wave function. We begin by rewriting 
equation (pd|) (in units of ~^^T) in a slightly different form 



where we have set 



¥(a;) = {er7 + »73^ + &(a0}*(a0 (29) 



7 = 7o cosh 9 — 74 sinh ( 



£7 = ; — e r cosh 6* 
V2A 



Pt = — e r sinhfl, (30) 
V 2A 
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1, 2 and e r is given by 

e, r = 1 

(31) 

-e, r = 2 

with e given by equation (|T7|), The parameter 9 is related to the particle velocity 
in the direction parallel to the wall, vt by vt = tanh#. 

From now on, we set pt = that is e = E (in units of -JL=T) so that the fermion 



moves along the direction normal to the domain wall. From ( p0| ) we also see that 7 
becomes 70. The general case for which the fermion has a non vanishing pt can be 
obtained by Lorentz boosting the solution along the direction of the wall by means 
of the transformation 

S(9) = exp(- 7Q7 t0). (32) 



To finish this section, let us work with the e-independent part of equation (|2 
and evaluate 

{+273-^- + tg(x)}$(x) (33) 
ax 

using $ given by ( |l8| ) and (^), with given by (15), the above equation becomes 

E ^± u ± (34) 

with 

V4 = (Ai)^^ + (A s ± )^^ a) (35) 
where we have defined, using cf>± given by ( p5| ) and (p6|), 

V± = {=F^ + ^)}0±(a:). (36) 
To construct the full wavefunction we set e = -jzt-E and 7 = 70 in (|29| ) and 



/2A 

boost then the so found solution in the direction along the domain wall. We need 
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to know therefore how 70 operates on the spinors u. In our representation 



To 



/ 
-/ 



and thus 

70 u ± = u T . (37) 



Recall that each of the four terms in equation ( |lq ) is by itself a solution to 
equation (|l6|) and so for a given sign of the energy, equation (p9j), with <J?(x) replaced 
by any of its four components, represents four solutions from which only two should 
be independent. In appendix A we show that ip± are related to cftzp by 

<p { ± a) = 2(^±a)^ Q) 

= 2(^T«)^ +Q) . (38) 

We also notice that (j)±~ a \x) can be obtained from (x) by exchanging a with 



—a. From equation ( f29| ) together with ( pq ) , ( p7\ ) and (pq), we get the four inde- 
pendent solutions to the Dirac equation labeled by the indices s and r 

= (A s ) ( -- a \e r <p { - a) u s _+2(d + a) ( f> { : a) u s + ) 

+ (^ s )( +a )(e r <AV +Q) ui+2(e-a)^ +Q) (39) 

where the constants (A s )(~ a ^ and (yl s )( +a ) replace (A s )| ) r a ' ) and (A 5 )^"^ respec- 
tively and have to be fixed by normalization and the choice of initial conditions. 
Again the general case for which the fermion has a non vanishing pt can be obtained 
by Lorentz boosting ( |39| ) along the direction of the wall according to (^) . Starting 
from the solution (|39] ) we will discuss in the next section the scattering states of 
fermions traversing the domain wall. 
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3.2 The scattering states 

At the electroweak phase transition the temperature is larger than the typical par- 
ticle masses such as that of the Z boson or light mass quarks and although it can 
be of the order of the top quark mass, most of the fermions will pass by the domain 
wall with energies above its height. To understand the interaction of fermions of 
such energies with the wall we proceed to study the scattering states of the fermion 
wave function and to derive reflection and transmission coefficients. 

The two physical processes of interest here are either a fermion leaving the 
bubble or a fermion falling into the bubble. We call the solutions corresponding 
to the boundary conditions of these processes the type I and type II solutions 
respectively. Both processes are depicted in Figure 3. 

For the type I solution we have an incident fermion from the right (x — > +00) 
with an energy parameter e which is always larger than the height of the wall ( 
e 2 > 4£ 2 ). At the wall this fermion is scattered into an reflected wave going to the 
right and a transmitted wave going to the left (x — > —00). Overall the fermion is 
therefore represented by an incoming and reflected wave to the right of the wall and 
by a transmitted wave to the left. 

If a fermion is incident from the left (type II solution) it can approach the 
wall with energies above its height (e 2 > 4£ 2 ) or below its height (e 2 < 4£ 2 ). In 
any case there will be again the reflected wave running now to the left and the 
transmitted wave running to the right. But for energies below the barrier (e < 4£ 2 ) 
the transmitted wave has to tunnel through the wall leading to a decaying instead 
of an oscillating wave. 

In this section we start out by constructing the wavefunction corresponding to 
the boundary conditions of the type II solution and determine then the transmission 
and reflection coefficients for this case. We show then that we need to rewrite the 
hypergeometric functions as expansions around z = 1 to obtain the wavefunction for 
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the boundary conditions of type /. It turns out that the transmission and reflection 
coefficients for this case are identical to the ones obtained for type II. 

For a given pair of indices s and r, equation ( |39|) consists of two terms with dif- 



ferent asymptotic behaviors as x — > +00. In this limit, according to equation (27), 
the first term proportional to °^ behaves like exp(— 2ax) whereas the second 
term proportional to (f)^ 01 ^ behaves like exp(2ax). The boundary conditions ap- 
propriate to the description of particles crossing the wall from the symmetric to 
the asymmetric phase are as follows: At x = —00, corresponding to z = 1 (out- 
side the bubble), we require the solutions ( |39|) to describe two plane waves, one 
moving towards (incoming wave fy mc ) and the other away from the wall (reflected 
wave ty re f). At x = +00, corresponding to z = (inside the bubble), we impose 
that for energies such that e 2 > 4£ 2 there is only a single plane wave moving away 
from the wall (transmitted wave \]/ irans ) 5 while for e 2 < 4£ 2 the solution dies out 
exponentially. 

The above conditions require that (A s )(~ a ^ = in (H) and the wave function 
for case II looks therefore like 

(%) n = A n (eJ+ a) u s _ + 2(£ - a)^ a) u s + ). (40) 

The normalization constant Ajj is determined in the next section and turns out to 
be independent of the quantum numbers s and r. 

To compute the transmission and reflection coefficients, we need first to look at 
the behavior of the solutions as x — > ±00. 

If we take the limit x — > +00 of equation (^) we find 

(Kfir = A n (e T u s _ + 2(£ - a)u%)e 2ax (41) 

For x — > —00, we have first to evaluate the hyper geometric functions in the 
second equation of (^) at z = 1. Since these functions are defined as expansions 
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around z = they are ill defined at z = 1 and we need to use the identity [E 

t? i u \ T(c)T(c - a-b) 

2 Fi(a,b;c;z) = — — 2 Fi(a, b,a + b - c + 1;1 - z) 

T(c — a)T(c — b) 



c _ a _ b T(c)T(a + b-c) 



r(o)r(6) 

x 2-Pi(c - a,c - b,c - a - b + 1; 1 - z). (42) 



therefore, using equation. (|42| ) in ( |2q ) and considering the limit x — > — oo we get 
from ( flO|) after some algebra 

-z^- + 

™ _ A 7J r(l - 2a)r(-2/3) / e r uj 2K-q)u^^ 



1 r,H T(-a-P + C)T(-a-p-0\-a-p-C -a-(3 + d 

ref A II T(l-2a)T(2P) f e r uj 2(g-q)< \ 2gx 

1 rJ// r(-a + /3 + 0r(-« + /?-0 V-a + /?-e -a + /3 + d 

(43) 

We can check that the reflected wave can be obtained from the incident by exchang- 
ing j3 with —(3. 

To compute the reflection and transmission coefficients, we have to compute 
the ratio of the reflected and transmitted fluxes to the incoming one. It is thus 
sufficient to calculate the ratios of the normal components of the corresponding 
vector currents. 

The normal component of the currents associated with the plane waves that we 
have found are 

j 3 = ^73^ (44) 

where is any of ty trans ; \$ mc or \j/ re /. Prom our solution, equation (|40|), the normal 
components of the incident, reflected and transmitted currents are for type II 

r(l + 2a)r(2/3) 



(j}f) 3 = A\A n \ 2 ee r 



r(i + a + p + f)r(a + - 
15 



T(l + 2a)T(2(3) 



r(l + a- (3 + QT{a- P-Q 



UT/)3 = -MAn\ 2 ee r 
{jfrh = 8|^| 2 e r |a| 

(45) 

The reflection and transmission coefficients R and T are now the ratios of 
the reflected and transmitted normal currents to the incident one, respectively, 
projected along a unit vector normal to the domain wall. Using equations (pSf) and 



the identity [22| 



T(x)T(-x) = . } , (46) 
x sin to 



we find after some simplification 

— sin 2ira sin 2-7T/3 



T 



sin7r(£ — a — P) sin7r(£ + a + 0) 



R = sirnr^ + a-f3) sin7r(g - a + 0) 
sin7r(^ — a — ft) sin7r(£ + a + (3) 

From the above equations we see that for states with energy parameter e 2 = 4£ 2 , 
the transmission and reflection coefficients become and 1. For states with e 2 < 4£ 2 , 



a becomes negative and real. From equation (41) we see that j trans is identically 
zero and therefore the transmission coefficient T is also zero. Correspondingly, the 
reflection coefficient R for fermions with e 2 < 4£ 2 is one. 

Both reflection and transmission coefficients are the same for positive and neg- 
ative energy solutions and thus for fermions and antifermions. They are depicted in 
figure 4. In this and the following figures we choose two representative values for the 
parameter £. The heavier fermions are represented by £ = 4 which corresponds to a 
ratio of fermion to Higgs mass of 2. It is expected that the mass of the top quark is 
lying in this mass region, so that we can consider the results for £ = 4 to represent 
a top interacting with the domain wall. The other fermions, on the other hand are 
much lighter than the Higgs. The corresponding value of £ will be therefore very 
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small. We choose a value of £ = 0.6 to represent these particles which corresponds 
to a ratio of fermion to Higgs mass of 0.3. 

We plot in figure 4 the transmission and reflection coefficients as a function of 
a reduced energy parameter y = The top of the barrier is then at y = 1. We 
see that the interactions between the heavy fermions and the wall quickly die out 
if we increase the energy parameter. The light fermions on the other hand are still 
feeling the presence of the wall high above the top of barrier. 

The result for the reflection and transmission coefficient is unchanged if we 
boost the solution along the direction of the wall to go to non-zero pt- Since the 
boost-operator defined in ( |32| ) commutes with 73 

[73,S(0)] = O 



we see that the currents defined in (H4j) is invariant under S. This shows that 
also (|47|) is unchanged. 

We proceed now to the solution of type /. Here the transmitted wave has to 
behave like e~ 2 ^ x for x — * —00. But the solution (|39|) is not regular in this limit - the 



hypergeometric functions diverge. To resolve this problem we have to rewrite (39) 
in terms of hypergeometric functions defined as expansions around z = 1. 

There are two possible ways to obtain this modification of (|39|), One possibility 
is to restart from (|24[ ) and write the two independent solutions of the hypergeo- 
metric equation as hypergeometric functions defined as expansions around z = 1. 
Equation ( j24j ) becomes then (|25f|) 

X± P \z) = 2 Fi{a T ,b T ,a T + b mp - c + 1; 1 - z) 

x ( ±®{z) = {l-zf~ a T-H 2Fl (c-a T ,c -b T ,c-a T -b T + l;l - z) (48) 



and ( |26| ) is now 



= z a {l-zf 2 Fi{a^,a^ + b^-c+l;l-z) 

z a (l-z)- 13 2 Fi(c -a T ,c-b T ,c-a T -b T + l;l-z). (49) 



17 



Finally we repeat all steps leading to (|39|). The superscripts in the equations above 
indicate the behavior of these functions at x — > — oo corresponding to 1 — z = 0. 



We can deduce this behavior by comparing (4£) with (]2 



A second approach for rewriting (|39|) is to start out with the four independent 



solutions we found in (39), use then ( Eg ) and the identity (p5|) 

2 Fi (a, 6, c; 1 - z) = z c " a ~ b 2 -Pi (c - a, c - 6, c; 1 - z) (50) 

to rewrite the four solutions in terms of hyper geometric functions of the variable 
(1 — z) instead of z. 

Both approaches yield as result 

% = (A s )^\e r <p i ^ ) u s _-2(3^ /3) u s + ) 

+ (A s )(-^(e r ^ ) ui + 2/3^ ) (51) 

To impose the boundary conditions for the solution of type I we have to require that 
for x —* — oo only terms oscillating like e~ 2 ^ x survive. We thus set (A s )^ +/3 ^ = 0. 
The solution of type I is therefore 

(*»)/ = A/M^u! + 2/3<A ( T ,3) 0. (52) 

Continuing now analogously to equations ([II]) to (47) we find the transmitted wave 



(**)/ " =^°° (V s r )f ans = Aj( er ui + 2Pu s + )e~ 2l3x (53) 
and the incident and reflected waves for this case 

A I Y{l-2(3)T{-2a) ( e r ui 2/3u^ \ „ 



(54) 



V r " r(-a-)9 + 0r(-a-^-0 V-a-/3-C -a - /3 + £ ' 

e/ Aj£(l - 2/?)r(2a) / e r iij 2/?< \ 2tM 

1 r)l F(a-P + £)!> - /3-e)U-/3-C a-/3 + r' 
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Here we see that we can obtain the reflected wave from the incident wave by ex- 
changing a with —a. 

As before we can calculate the currents and find 

2 



\Ar\ 2 \a\e r 



T(l + 2/3)r(2o) 



U? f ) 



u 



trans \ 



-8\Aif\a\e r 
l/| 2 e r e 



r(i + a + /3 + <e)r(a + /?-£) 

r(l + 2/3)r(2a) 



T(l-a + P + Z)T(-a + p-Z) 



4\At> 2 - 



(55) 



If we compare these currents to the ones obtained in type II (45) we see that 
we obtain (|5^) from (p5j) by exchanging a with (3. Since the transmission and 
reflection coefficients in ( f47|) are unchanged under exchange of a with [5 we find that 
the reflection and transmission coefficients for a fermion falling into the asymmetric 
phase or a fermion escaping the asymmetric phase are the same for energies above 
the barrier and are given by (|47|). 



4 Normalization and Orthogonality 



In this section we want to construct a general solution to the Dirac equation using 



the type / solution given in (52) and the type II solution given in (f40[). The general 
solution should have the form 



^ S r = a s r (^ s r )i + K(n)n 



were o, b are arbitrary coefficients normalized to 1 



l = (a s r ) 2 + (Kf 



(56) 



(57) 



In (|56| ) and \£ jj should represent then the normal modes of the system as depicted 
in figure 3. 
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To obtain such a solution we have to normalize first and themselves and 
then assure that the two solutions form an orthogonal set. It will turn out that they 
actually are not orthogonal and we finally will show how to orthogonalize them to 
obtain a solution of the form ( |56| ) . 

To normalize the type I and II solutions given in (|52| ) and in ( |40| ) we first choose 
an open interval (+/, —I) on the x-axis of figure 1. We define now the normalization 
condition for a state on this interval as follows 



limi^oo I dx^^. (58) 



For I->oowe recover the usual normalization condition from (| 

To evaluate the integral in equation (^) we introduce an additional length scale 
5 on (— I, +1) such that 5 <C I. The scale 5 defines an symmetric interval around 
x = such that inside [—5, +5] the Higgs potential g(x) changes rapidly from its 
value close to zero, to its value close to 2. 

rl r-S r+S r+l 



f dx^y = ([ + f + I ) dx^m 
J-l J-l J-S J+8 



If the scale I becomes very large we can replace the wavefunction ^ on the intervals 
(— I, —5) and (+<5, +1) with its asymptotic limits at x — > =Foo respectively and obtain 
so 

£ ^* f * = J t dx^^^ + J +l dx*\ +oo) * {+oo) + 0(5/1). (59) 

Here ^±oo represents the asymptotic limit of the wavefunction at x — > ±oo respec- 
tively and the integrals have to be evaluated only at the indicated limits ±Z. If we 
let I — ► oo the corrections of order 0(5/ 7) vanishes and expression ( |59| ) becomes 
exact. The result diverges in this limit linearly in I analogous to the divergence in 
the volume of the plane wave normalization. 
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To solve equation ( |59| ) for the constants Aj and Ajj of and ^77 we use their 
asymptotic expressions given by (|53|), (pH), fl4l| ) and (p3| ) respectively. Using the 
orthogonality of the spinors u in (|l5|) and the fact that integrals over oscillating 
functions average to zero / dxe lkx = we obtain after some algebra 



4k 2 \ + r + Ht 



a ' 



As before we can obtain Aj from Ajj by exchanging a with f3. 

To prove the orthogonality of ^/ and Vl/jj we need to evaluate their overlap 
integral given by the scalar product 

1 = J + dxV\Vn (61) 

The integral is defined as before on an interval (—/,+/). If the overlap integral / 
in ([H]) is zero then and ^S>n are orthogonal, if it is non-zero we still have to 
search for orthogonal solutions. 



Proceeding analogously to equations ( 5q) to (|59|) we obtain 

1 = j + da;(*J) (+oo) (*jj) (+oo) +^ ^(*l)(-oo)(*//)(-oo) (62) 
Using again the expressions for the asymptotic wave functions we find 

r s/4 4Vw n , r(2/3)r(-2a) 

r(/3 - a + £)(/? - a - 4 + 1) 
# (63) 

The overlap integral is therefore non-zero, ^1 and ^ n are not orthogonal. To 
obtain a solution containing normal modes like it was envisioned in (|5^) requires 
therefore to find a set of orthogonal states. This can be done using the Schmidt 
Orthogonalization (pGfl). 

In the Schmidt Orthogonalization we start out with a set of non-orthogonal 
states and construct out of these one by one an orthonormal basis. The set of 
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non-orthogonal states is in our case Vl/j and Vl/jj. Let's choose now tyj as our first 
normalized basis state. The second orthonormal basis state corresponding to tyj is 
then obtained through the ansatz 

Martha = N(aV i + V n ). (64) 

Here N and a are constants to be fixed by the requirement that ^> rtho is normalized 
to one and is orthogonal to Vl/j respectively. 

The orthogonality condition is evaluated by forming the scalar product of ^ortho 
in (H) and *j 

J dx^ortho = N(a J dzttjtfj + J dx^n) 

= (65) 

The first integral on the right hand side is one since \Pj is normalized and we obtain 
therefore 

a = -fdxV\9 u = -I. (66) 



The constant a is therefore identical to the overlapp integral / evaluated in (|63|). It 
is now straightforward to normalize ^ortho- The normalization condition reads 

1 = J dx*i tho * ortho (67) 

which can be solved for 



We therefore found a general solution of the form ( p6[ ) with Vl/jj replaced by ^ or tho- 
Of course, we could have started in (|64|) with the orthonormal basis state to 

9' ortho = N'(a'* II + *i). (69) 

It is straightforward to see that in this case a' = — I* and \N'\ 2 = \N\ 2 , so that up 
to a phase ^' ort h = ^ortho- Again we can construct the general solution of normal 
modes, this time replacing ^/ in ( |56| ) with ^* rtho - 
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5 Conclusion 

In this paper we discussed the motion of fermions under the influence of electroweak 
domain walls. We showed that such fermions are well described by a Dirac equation 
with an effective mass term. This mass term was obtained via the Yukawa coupling 
in the Lagrangian applying the classical mean field approximation. We substituted 
for the finite temperature vev of the Higgs field the classical solution to the equations 
of motion of a Higgs field in a finite temperature effective potential. 

We investigated the Dirac equation analytically and found the wave function of 
the fermion. Transmission and reflection coefficients were derived and found to be 
the same for both fermions and antifermions. 
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Appendix A 



Here, we would like to show that 92's defined as 

tp± = [=F^ + £s(zO]0± (i) 

in equation (|36|) are proportional to the (f> 's given by equation (p6|) where 5(3;) = 
1 + tanhx. We will prove it for one case, show that it can be done for the others 
and write the specific relations for all cases. 

Let's start by looking at one case, for example, 

= [-^ + &(*)]*t a) (ii) 

where 

^+ a) = z a (l- zf 2 F 1 (a + (3 -f + l,a + /3 + e,2a + M) (iii) 
and z = |(1 — tanhx). Then it follows that 

[-^+ig{x)]=2(l-z)[z±+£] (iv) 
so we can rewrite tp{ a ' as 

<p$ a) =2{l-z)[z±+t]<fc a \ (v) 



Let's consider 



Multiplying both sides of this equation by 2(1 — z)z we obtain 
2(1 - z)z~^ 1 ^ ^ a) ] = 2(1 - Z )[z± + ^]^ Q) - 



(vi) 
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Comparing this with equation (26) and using the explicit form for cj), we have 



</4~ Q) = 2(1- z)z-^ +l ^ +a {l- zf 2 Fi (a + /?-£ + 1, a + /? + £, 2a + 1,*) 

az 



= 2(1 -z)z a+1 [—(l-zf] 2 F 1 (a + P-S + l,a + + S,2a + l,z) 
+ 2(1- zf +1 z-t +1 -^-[zt+ a ^(a + /?-£ + l, a + /3 + £, 2a + M)] 

We now add and subtract (3 to the exponent of z inside the bracket in the second 
term above and differentiate to get 

f z [z aH 2 Fi(a + (3-i + l,a + l3 + i,2a + l,z)} = 
z~P £ [ z a +P+t 2Fl (a + (3-Z + l,a + (3 + Z,2a + l,z)} 
- (3z a+ ^ 1 2 Fi(a + P-Z + l,a + P + Z,2a + l, z). 

Now we use the identity ( J25| ) 



-^[z b 2 F 1 (a,b,c,z)\ = bz h ~ x 2 Fx(a,b + l,c,z) 

and putting everything together, we find 

(f { + a) = -2[3z a (l- zf 2 Fx(a + P-^ + l,a + P + ^,2a + l,z) 

+ 2(a + (3 + t)z a (l- zf +l 2 Fi(a + /3-£ + l,a + /? + £ + l,2a + l,z) 

In the second term of this equation we can use the identity (p5|) 

c — ct c — b 

2 Fi(a,b,c,z) = — — 2 Fi(a- l,b,c,z) - — — 2 Fi(a,b- l,c,z) 

(b — a)(l — z) (b — a)(l — z) 

and then combine the terms such that 

<P+ a) = [-W--M 2 -(l3 + if}]z a (l-zf 2 F 1 (a + + l,a + P + £, 2a + l,z) 

+ \[{a + i) 2 -P 2 ]z a (l-zf ^(a + P-^a + P + i + l^a + ^z) 

Using the fact that a 2 = /? 2 + £ 2 , we see that the coefficient of the first term is 
identically zero and that 

j[(a + t?-f3 2 ] = 2(t + a) 
25 



Finally comparing to equation ( |2"6| ) 
we see that 

^~ a) = 2(e + a)^ ( : Q) (vii) 

Similarly we can show that 

^ +Q) = 2(e + a)^ +Q) (viii) 

For the case when yj's and 0's are expanded around (1 — z) instead of z, we can 
apply the same method to derive similar relationships. They are 

^ = -2/3<t£® (ix) 

These relations are used in the calculation of the fermionic wave function to elimi- 
nate the linearly dependent solution. 
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Figure captions 

Figure 1: The bubble profile, x = — oo corresponds to the bubble's outside. The 
domain wall is the region where the vev of the Higgs field g(x) changes rapidly. 
Figure 2: The dimensionless effective potential for the Higgs field V(g, Q in units 
of 5T(^) 3 for different values of ( = 0, 2, 1.75, 2.25, 2.75. 

Figure 3: Scetch of the asymptotic behavior of the solutions of type / and 77. 
Figure 4: Transmission and reflection coefficients for £ = 0.6,4 which correspond 
to a ratio of fermion to Higgs mass of 0.3 and 2 respectively. The reduced energy 
parameter is y = so that the top of the barrier is at y = 1. 
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